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Abstract 

The leading correction to the smoothed connected energy density-density 
correlation function is obtained for the large energy difference, within the con- 
text of the Gaussian Random Matrix Theory. In order to achieve this result, 
the supersymmetry method is extended beyond the sigma-model, to include 
small quadratic fluctuations around the saddle point. Special care is taken to 
avoid the potential divergence arising from the unbounded nature of the saddle 
point. Also, in the small energy difference regime, the leading correction to a 
two point correlation function is obtained. 

Keywords: Random matrix theory, Supersymmetry, correction. 
PACS-96: 05.40.+j 



1 



1 Introduction 



The Random Matrix Theory (RMT) has its roots in an idea of Wigner that statistical 
mechanics should be used to analyze highly excited energy levels of large nuclei 
This concept was further developed by Dyson p[ to become what is now known 
as RMT. Since then, it has been applied to numerous problems related to complex 
quantum systems in different areas of physics f^. 

In RMT, the Hamiltonian H is represented by a matrix of dimension A^, for N 
large. The matrix is completely random, except that it satisfies a certain symme- 
try. The ensemble average of any function f{H) of the Hamiltonian is calculated by 
multiplying the function by a probability distribution and averaging over all possible 
matrices that possess the required symmetry: 



RMT has been utilized for various purposes. For example, in condensed matter 
physics of disordered systems, the Hamiltonian, or a part of it, can often be taken to 
be random, except for some physical symmetries. Assuming that a particular physical 
quantity does not depend on the details of the system, it can be calculated using RMT. 
The problem of dealing with a complicated (or unknown) Hamiltonian of the physical 
system is replaced in RMT by averaging over an ensemble of random matrices which 
possess a certain symmetry, as required by the symmetry of the original Hamiltonian. 
Commonly used symmetries of the ensemble are unitary, orthogonal and symplectic. 
The problems where RMT has been applied to include conductance fluctuations in 
mesoscopic systems persistent current of disordered metal rings |^ and the motion 
of a quantum particle in a thick wire [0. A similar approach is taken in analyzing 
the low energy behavior of QCD, only here the ensemble considered possesses one of 
the chiral symmetries: chiral orthogonal, chiral unitary, or chiral symplectic 0. 

2 



{f{H)) = -j VHf{H) exp{-ArTr [V{H)]} 



(1) 



where Z is the partition function: 




(2) 



Another incentive for exploring RMT comes from its relation to both classical and 
quantum chaotic systems. On the basis of numerical evidence, it was strongly conjec- 
tured that certain statistical properties of quantum chaotic spectra are universal and 
identical to those predicted by Gaussian ensembles of RMT The connection was 
further elaborated in and proved in semiclassical limit in |10|. Another important 
problem in chaos, the transition from classical integrable to chaotic behavior, has 
been investigated by relating classical chaos to the size of a band of an ensemble of 
banded random matrices |Tl|. Banded ensembles have been constructed to describe 
the level fluctuations of a particular type of anharmonic oscillator with a parameter 



that is responsible for the transition from integrable to non-integrable behavior 

In all the applications of RMT mentioned above, one deals with Gaussian ensem- 
bles. The potential V in (|I|) has a simple quadratic form (^), which is equivalent 
to the requirement that different elements of matrix H be uncorrelated. It is in- 
teresting to investigate how additional terms in the potential of the non-Gaussian 
ensembles affect different correlators, or, in other words, what are the quantities cal- 
culated within RMT that do not depend on the form of the potential, but on the 
symmetries alone. Universal properties of correlations for non-Gaussian ensembles 
and their relation to the Gaussian ensembles were studied and the correlations calcu- 



lated, in the context of matrix models, in [13| using the loop equations method. The 



same properties were investigated in the context of disordered systems in |T^ using 
the orthogonal polynomials. The conjecture made there was clarified in [|T^, while 
the rigorous treatment was given in [l^. Other methods, such as the diagrammatic 
method |l7|, the functional derivative method |jl8|| and the supersymmetry method 
[ P^ were also introduced and used for similar problems. 

On the other hand, RMT can be extended to matrix models which have been 
studied for their own sake. Matrix models, for example, has been used to describe 
2D gravity. The partition function of the matrix models is obtained from the RMT 
partition function by adding a kinetic term to the potential and summing over all 
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Feynman diagrams, while in 2D gravity, the partition function is obtained by sum- 
ming over all surfaces. The potential term generates Feynman diagrams, which can 
be regarded as the inverse lattice of a discretized surface of arbitrary genus. More 
precisely, a Feynman diagram is the inverse lattice of the discretized surface. The 
ensemble average of RMT can now be regarded as a sum over all diagrams (surfaces). 
Furthermore, when is large, the partition function has a perturbative expansion in 
A^~^, each term corresponding to the contribution of all surfaces of a fixed genus. It 
has been shown that the dominant contribution comes from a sphere (genus 0) ||20|| , 
so the approximation in which only the leading term in the expansion is used is called 
spherical or planar approximation. The genus 1 contribution to the partition function 



and the one point correlator have been calculated in using the loop equations. 
Good reviews of different aspects of matrix models are given in [ P^] . 

While the original approach used in RMT was the method of orthogonal polyno- 



mials 0, the supersymmetry method, pioneered by Efetov |2^, has been also proven 
to be very useful, especially when the coupling to external sources is needed M . It is 



essentially an extension of Anderson's replica trick p4|. A very thorough and detailed 
review of the method can be found in |25[ . 

The results cited above were concerned with the leading order contributions to dif- 
ferent RMT correlators. The leading order correction (genus one contribution, in the 
language of matrix models) to the one-point correlation function for the Hermitian 
matrix model has been calculated in PB[, while the iterative procedures to calcu- 



late other corrections (higher geni contributions) were presented for Hermitian ||27 



complex 1 2^, and orthogonal and symplectic |^ matrix models. Recently, there has 
been significant interest in finding the corrections within the context of disordered 
systems. Orthogonal polynomials have been used to calculate leading corrections 
to one and two-point correlation functions for the complex RMT The super- 

symmetry method has also been used to calculate corrections to various quantities 
within RMT: the energy-energy correlation in the regime where the energy difference 
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is small [0, the level curvature distribution and the eigenfunction amplitudes 
distribution p^ . 

In this paper, the leading order correction to two-point correlators, in the A^~^ ex- 
pansion, is calculated using the extension of the supersymmetry method. Within the 
usual supersymmetry method, the partition function is expanded around the saddle 
point (the saddle point manifold, to be precise), the energy difference is considered 
small (~ N^^), and only the contribution coming from the saddle point is considered. 
This simplification leads to the so-called zero- dimensional non-linear a-model. Here, 
the method is extended beyond the cr-model: both to accommodate for the large en- 
ergy difference (~ N^) and to include small fluctuations around the saddle point. By 
extending the method for the large energy difference the leading contribution to the 
two point correlators can be obtained (in the large energy regime), while the inclusion 
of the fluctuations allows the calculation of the corrections in both regimes. As it 
is customary in the large energy regime, the smoothed energy density-density corre- 
lation function is considered {Ti]; that is, the correlator is averaged over the energy 
scale of the order of few A^~^'s, so that the terms oscillating with the frequency of 
the order go to zero. The leading order contribution of the (smoothed) connected 
energy density-density correlation function obtained in this text agrees with the result 



obtained using the method of orthogonal polynomials, [|IJ]. In the small energy differ- 
ence regime, the correction to the connected two point correlator, {G"^ {Ei)G~ {£2))'^ , 
is obtained. The leading term of the correlator is also calculated and it is found to 
agree with the well known result which was obtained using the supersymmetry, as 
well as other methods. 

Furthermore, the convergence of the method is proved. Particular care has to be 
taken because of the potential divergence arising from the unbounded nature of the 
saddle point manifold. The singularity structure of Green's functions is used in order 
to resolve the ambiguities arising from the double- valued nature of the saddle point. 
As it is usual in the supersymmetry method, the Gaussian ensemble is considered in 
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the limit when N ^ oo and the energies are taken not too close to the ends of the 
spectrum. For the sake of the definiteness, in this text, the ensemble with the unitary 
symmetry is considered. 

The approach used in this text was to follow the presentation of the supersym- 
metry method as given in an excellent introductory paper by Zuk [3^] and to modify 



it when needed. Some parts of the calculations were done using the symbolic soft- 
ware package Dill as well as additional programs written in the Mathematica© 



programming language [36 



The paper is organized as follows: the important steps of the usual supersymmetry 
formalism are given in section |^. The definitions and conventions related to the 
Grassman variables and the integration over the Grassman variables are given in 
appendices ^ and ^, respectively. The formalism is modified to include the corrections 
in section while the correlation functions are calculated and the convergence of the 
method is proved in section |, using the results presented in appendix |C|. 

2 Short Review of the Supersymmetry Method 

In this section, the main points of the supersymmetry method are summarized. This 
is needed in order to show how the extended supersymmetry method is constructed 
and to prove its convergence. We begin by defining the relevant quantities. The 
quantity of main interest here, the connected two point correlation function, is given 
by: 

(G+(Ei)G±(E2))^ = {G+iE,)G^iE2)) - (G+(Ei))(G±(E2)), (3) 

where: 

G^(E) = Tr , (4) 

and the averaging (. . .) is performed over Gaussian unitary ensemble; therefore, it is 
given by (|l]) with the potential V given by: 

n^) = ^^^ (5) 



where A is an arbitrary constant. 
The density of states is given by: 



P 



2tiN 



G+{E)~G-{E) 



(6) 



so that the connected part of the ensemble averaged energy density-density correlation 
function is given by: 



{p{E,) p{E,)) - {p{E,)) {p{E,)) 
1 



27r2Ar2 



Re 



{G+iE,)G-{E2)r - {G^{E,)G^{E2)y 



(7) 



The density of states satisfies Wigner's semi-circle law 

1 



p{E) 



2'kX 



y4A2 - E\ 



(8) 



so that the energy eigenstates of H are confined to a region between 2A and — 2A. 

In order to begin with the supersymmetry formalism, let us consider now a gen- 
erating function: 



= exp{— Trg In [D + J{s)]} 



(9) 



where Detg and Trg denote the graded (super) determinant and the graded (super) 
trace respectively and are taken over all indices. The symbol trg is used when the 
summation is to be taken over p and a indices only. (For a complete introduction to 
anti-commutative (Grassman) variables one can consult and |^.) The definitions 
and conventions used in this article related to the Grassman variables are given in 
appendix D and J are 4N x 4N graded matrices, diagonal in both upper and 
lower indices, the only non-vanishing elements being: 



D 



00 

pp 



11 
pp 



{El + iep) In - H, 
(-)%ljv. 



(10) 
(11) 
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It is then easy to show that: 



G{Er+ie^)G{E2 + ie2) 



e=0 



-Zie] 



(12) 
(13) 



£ = 



4 deide2 

The two point Green's functions, G^{Ei)G^{E2) and G^ {Ei)G^ {E2) , are obtained 
from the above equation by taking ei = e = —£2 and ei = e = 62 respectively, for 
e positive and infinitesimaL Therefore Z{e) is indeed the generating function for 
the two point Green's functions. In this case, the matrix D + J consists of two 
diagonal blocks, corresponding to the two Green's functions. The two blocks contain 
the information about the respective energies as well as the sides of the branch cut of 
the Green's function from which the energy approaches the real axis. If the ensemble 
averaged product of more Green's functions were to be calculated, additional blocks 
would have to be added to D + J. 

Some important properties of the integration over Grassman variables are given 
in appendix |B[ Using the Gaussian super- integral ([73|) , the generating function (^, 
can be ensemble averaged and shown to be: 



Z{e)) = / V(PV(P (exp 



p,a,fj, 



exp I 



^ 0^(/i)(E + ^L + ie,h + ie2h + JWM 



(14) 



where the average energy E = {Ei + E2)/2, and the difference between the energies 
uj = El — E2 have been introduced and L, Ii and I2 are diagonal in p and a indices: 



(15) 



At this point, it is clear that due to the Grassman variables, the generating func- 
tion is normalized by the definition (^ (^(0) = 1) which makes it possible to in- 
terchange the integration over 0's and the ensemble averaging (integration over the 
random matrix H) and, eventually, to perform the ensemble averaging completely. 
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In the case where a generating function is not normahzed, it has to be normahzed 
first by dividing by the corresponding partition function and then ensemble averaged. 
Typically, both the partition function and the generating function are integrals, so it 
is a very nontrivial task to perform the ensemble averaging. 
The procedure is continued using: 



(exp 



X 



) = exp -—trg 5^ s^^: = y: r.Mi'if^). m 



Then the ensemble average over the random Hamiltonian can be completely performed 
to give: 

A2 



{Z{e)) = J l?00exp 



2^-trg + itrg {E + ^L + ieJi + ^£2/2 + J)S 



;i7) 



Obviously, the price that has to be paid for carrying out the ensemble averaging 
completely is the introduction of the quartic terms (S"^) in the integrand. The way to 
proceed is to eliminate these terms, using the Hubbard- Stratonovich transformation, 
as follows. 

The function W is defined by: 



iN 

+ trg (ei/i + e2/2)(^ - —a) 



(18) 



where a is an arbitrary 4x4 super-matrix having the same symmetry properties as 
5*. Using the parametrization for a introduced below, and as long as is large, it 
follows that: 

J Pae*^('^'^) = J I)ae*^('^'°) = 1. (19) 

The second identity is easy to check once the parametrization for a is introduced 
(P^-(P^). After the previous identity is inserted in (|T^ and the shift of variables: 
(7 a — j{^L + J) is performed, the averaged generating function takes the form: 



\Z{e)) = J Da exp |-ytrg a - —trg aL - —trg [eih + e2l2)(r- 



2A 

A^trg \n{E - Xa) - —trg 
A 



A 



a + + jieJi + €212) 



(20) 



J 



completing the Hubbard- Stratonovich transformation. 

At this point, it is important to notice that not only have the quartic terms dis- 
appeared from the generating functional, but in addition the number of integration 



variables have decreased from AN complex variables in (|Tj) to those spanned by the 
4x4 complex super-matrix a in (pOl). Not all entries in a are independent, because 
it has to satisfy certain symmetry properties. Therefore, in order to see what exactly 
the integration variables are, the symmetries of the super-matrix a have to be investi- 
gated first and then the parametrization of a (the parameters being the independent 
variables) consistent with the symmetries can be found. The detailed account of 
the symmetries and the parametrization leading to the independent variables of the 



super- matrix a for the Gaussian orthogonal ensemble is given in p5|. The analysis for 



the case needed here, the Gaussian unitary ensemble, is similar and is clearly stated 
3^ ; what follows is only a short description of the results. 



m 



In order to satisfy (0), and having in mind the definition of W given in (|T8D, the 
super-matrices a and S have to have the same symmetry. From the definition of S, 
([IBD it follows that: 5^ = sSs, so: cr^^ = sas also, leading to the decomposition: 

a = T-'PdT, (21) 

where 4x4 super-matrix T has to satisfy: 

= sT^s, (22) 

and Pd has to be a diagonal super-matrix with real boson-boson and imaginary 
fermion-fermion elements. The last condition has to be changed slightly in order 
to achieve convergence of integrals in (|T^) and make Pd belong to the saddle point of 
the exponent in (|20|) . Instead of (plD , a is decomposed as: 

a = T~\ao + N'^/^Pd)T. (23) 

The change amounts to a simple shift of the diagonal matrix Pd by a constant diagonal 
4x4 super-matrix ctq determined below. Pd now has real boson-boson and imaginary 
fermion-fermion elements. Also, Pd has been rescaled for future convenience. 
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Due to its defining property (^2]), the boson-boson part of tlie super-matrix T 
spans the (non-compact) group SU{1, 1) while the fermion-fermion part spans the 
(compact) group SU(2). Further, T can be written as T = RTq, where [R, L] = 0, so 
that R is block-diagonal in the lower indices, and the blocks are denoted by -Ri and 
i?2- The symmetry of the matrix T implies: T^^ = sT^s, R{^ = R[ and -R2 ^ = kR\k. 
This gives the final decomposition of cr: 

a = T^\cTo + N''/^6P)To = T-\PTo, (24) 

where 6P = R^^PdR- Thus 5P is a block-diagonal super-matrix with real boson- 
boson and imaginary fermion-fermion eigenvalues. The decomposition of a ([2^), can 
be regarded as splitting of the manifold spanned by the independent variables of a 
to the space spanned by 5P and the coset-manifold defined by Tq. Due to the change 
of variables, the original measure Va becomes: 

Va = I{P)VPVfi{To), (25) 

where Vfi{TQ) is the invariant measure on the coset manifold and I{P) is the Jacobian 
of the transformation I{P)'- 



(A?-A0)(A1-A^) 



1 2 



(26) 



_(A?-Ai)(Al-AO) 

and A"p are the eigenvalues of Pp. The derivation of the last expression is straight- 
forward; useful hints are given in similar derivations of and | |25[| . 



Following its symmetry, the matrix Tq can be decomposed as: 

To = U-'TM, (27) 

where: 
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T] and p are Grassman variables, and 












VI + A«oA«o 



(29) 



= -e^-^o sinh ^, //q = -e-^"^" sinh f 



Pi = ie"-^^ sin — , pi = ie 
2 



(30) 



6o, 01 and are bounded (from to 27r), while is unbounded. 
It is also useful to introduce: 

Ao = cosh 9q — 2popo + 1, 1 < Aq < oo 
Ai = cos 6*1 = 2pipi + 1, — 1 < Ai < 1. 



(31) 



With the change of variables stated above, the invariant measure on the coset manifold 
is found to be: 

1 



■Dp{To) 



- dXodXid(f)od(f>id'r)d'r)*dpdp* 



(32) 



(Ao - Ai)2 

It is obvious from the previous equations that the symmetries of the super-matrix 
Td are SU{1,1)/U{1) for the boson-boson part and SU{2)/U{1) for the fermion- 
fermion part. Some parameters of the original SU{1,1) x SU{2) of T have been 
absorbed in the degrees of freedom other than those of T^, nevertheless the non- 
compact degree of freedom can be identified as Aq. 

Finally, due to the symmetries of its factors, SP is decomposed to: 



where: 



exp 



SPi 
5P2 — exp 



SP = 



Ci* 
-Ci 

tQ 
-K2 



5Pi 

5P2 

ai 

ihi 

a2 

ib2 



(33) 



exp 



exp 



-cr 









Ci 

-^C2 

<2 



(34) 
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(i and (2 being Grassman variables and oi, 02, 61 and 62 ranging from —00 to 00. 
The measure on 6P is easily computed to be: 

-1 



VP 



;dai da2dbi db2d(id(* ^ d(2dC 2 ■ 



(ai - ihiY{a2 - ih2y 
Using the above stated parametrization for 5P and the definition of I{P) 



(35) 
it 



follows that: 



m 



((^01 - ^02 + N-^l^{ai - a2))(aoi - 0^2 + N-^l^i{hi - 62)) 



(36) 



.(^01 - (^02 + Ar-i/2(ai - ih2)){(Jm - ^02 + N~y^{ihi - a2))_ 

where the a^p are the saddle points determined in the next section. Combining (pSj), 
( |32D and (p5|), the final expression for the measure on the a-manifold is given by: 

dXodXid(j)od4>idr]dr]* dpdp* daida2dbidb2d(id(* id(2dC* 2 



Va = -I{P)- 



(Ao - Ai)2(ai - i6i)2(a2 - ^62)^ 



(37) 



with I{P) given above (36). 

Therefore, an exact integral expression for the ensemble averaged partition func- 
tion (^) has been found in this section. Also determined is the parametrization of 
the super- matrix a in terms of the independent integration variables (p^), (p8D , (p9D, 
(PID, ( P5D and (0), and the expression for the measure in terms of the integration 
variables (37). The connected part of the ensemble averaged density-density correla- 
tion function, the quantity of the main interest in this article, is determined from the 
derivatives of the averaged partition function, as follows from (pt) and (ITBp. 



3 Super symmetry Formalism Beyond the cr-model 

In this section, the calculation of the averaged partition function is continued using 
the saddle point method. Contrary to the usual supersymmetry method, where uj is 
taken to be small (of the order of A^~^) and only the highest order terms in are 
retained, both the leading and some sub-leading order terms are kept in this article. 
In other words, in addition to the contribution coming from the saddle point as in the 
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standard supersymmetry method, the small quadratic fluctuations, massive modes, 
are also considered here. 

The saddle point of the exponent in ( |20D can be obtained as a solution of the 
following equation: 



_5_ 

6^ 



itrg + l^trg aL + jtig {eJi + €2/2)0" + trg \n{E - Xa) 



0. (38) 



The solution of the previous equation is not unique. If the saddle point is required 
to be diagonal, the diagonal elements aop have to satisfy the following equation: 

The solution, for both values of p, is given by: 



aop = E' + i-Yu' - 2te'p ± i^l - [E' - {-)puj' + 2ie'p], (40) 
where E, u and has been rescaled: 

^Opi {p = 1, 2) can be regarded as the two factors of the correlator {G{Ei)G{E2)). 
Furthermore, for ep positive or negative in (^), G^{Ep) or G^{Ep) respectively, are 



obtained, as shown below (p^ . The Green's function G{z) as deflned by (|D {z 
complex), has a branch cut for z real and — 2A < z < 2 A. That means that the 
points G{E + ie) = G+{E) and G{E - le) = G-{E), for -2A < E < 2 A, belong to 
different Rieman surfaces, so that the sign of the small imaginary part ie determines 
the correct surface. On the other hand, the square in (^) also has a branch cut. 
Furthermore, it is again the sign of the e'p under the square root that determines the 
sign of the imaginary part of the square root and consequently the Rieman surface. 
Therefore, the branch cut structure of the Green's function arises from the branch 
cut of the (square root in the) expression for the saddle point (Xop- As in the standard 
cr-model, the convergence can only be achieved if the ± sign in front of the square 
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root corresponds to G^{Ep) p5|. The convergence of the extended model, presented 
here, is proved in the next section. Therefore, the solutions of ( pQ] ) are: 

a+T = E' + J-e' — ±iJl-iE' -uo'f±i^ (42) 

^l-{E'-uj'Y V V ; 

Eventually, e (in the equations above) will be set to zero, but the convergence prop- 
erties have to be determined before. Further, the diagonal saddle point is defined 
as: 

o-Q = diag(aoi, ctqi, 0-02, (Tq2)- (43) 

In the standard supersymmetry method, uo is taken to be small, so that the second 
term in the equation for the saddle point, (pH), is neglected. In that matrix of 

the form Tq^^cto^o also satisfies the saddle point condition. This is because the first 
and the last terms are obviously unchanged by the transformation oq ^ Tq^oqTq and 
the third is small because it is proportional to e. Consequently, the saddle point is 
actually a manifold, parametrized by the degrees of freedom of matrix Tq. 

The approach taken in this paper is that uo is not neglected. Using the parametriza- 
tion introduced in the previous section: (PB]), (^) and (^Tj), the transformation of the 
second term of the saddle point equation (|38|) can be calculated: 

trg Tq-VoToL = (fioi - ao2)(Ao - Ai) ^ trg a^L. (44) 

This lifts the degeneracy of the saddle point. If uj is small, uj ~ A^~^, the manifold 
spanned by the matrix a is almost fiat in the direction of the matrix Tq compared 
to the directions of the other degrees of freedom (matrix 5P). Consequently, the 
manifold Tq^^ctoTq can still be considered to be a saddle point manifold. Therefore 
it makes sense to integrate over it (over the almost flat direction of the matrix Tq), 
provided that the following constraint is observed. Since Aq, one of the degrees of 
freedom of the matrix Tq, is unbounded ([25|)-(PTD, its upper bound should be small 
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compared to in order that the second, degeneracy hfting, term of the saddle point 
equation (^) stays small compared to the other terms. The upper bound of Aq is 
still taken to be infinity, provided that in the limiting procedure N ^ oo first and 
only then the limit of the upper bond of Aq is taken. 

In the large u regime, the degeneracy of the saddle point manifold is lifted fully, 
so that the degrees of freedom of both Tq and 6P can be considered as fluctuations 
around the saddle point. 

The decomposition of the matrix a given in (p^ should now be much more trans- 
parent: it is simply the expansion around a saddle point. N^^^'^Tq^5PTq are the 
fluctuations around the saddle point and Tq represents the saddle point manifold or 
the fluctuations, depending on the regime. The factor N~^/'^ is introduced because it 
is customary in the saddle point method to suppress the fluctuations by the inverse 



square root of the coupling constant In effect, it gives another limiting procedure 
- the infinite bounds of the bosonic degrees of freedom of the matrix 5P are of the 
order N^^'^. 

Both (G'+(Ei)G+(^2)) and {G+ {E{)G- {E2)) are calculated from {Z{e))- the dif- 
ference lies in the sign of 62 and in the saddle points (^2]), used in the expression 



for {Z{e)). Neither ei nor 62 play any role in the further calculations related to 
{G^ {Ei)G^ {E2)) , so they can be set to zero. However, they are still needed for 
{G^ {Ei)G~ {E2)) . Both facts can be expressed using uj^' with the understanding 
that: 

{uo' + it' for the calculation of (G^ [EAG^ {E2)) 
(45) 
uj' for the calculation of {G+ {Ei)G+ {E2)) 

Expanding the exponent of the averaged partition function, (|2D|), around the sad- 
dle point (plD, keeping the terms to the order of A^'' in the expansion of the logarithm, 
and bearing in mind the definitions for (Tq (^31) and J (pA|), it follows that: 



{Z{e)) = J Vaexp{-2Nuj+'tTg T^^ao + N-^''^5P)TqL - 



\rg [1 - ((To + 2uj+'Lf] {5Pf + 2A^^/2^+'trg L5P 



2 
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^trg \T^\a^ + N-"HP)To + 2u^' l\j + 0{N-^'^)}. (46) 

The above equation makes it clear that the degrees of freedom of o", contained in 
the parametrization of 5P, correspond to the massive fluctuations around the saddle 
manifold. Only the quadratic fluctuations are kept in the previous equation. The 
others are suppressed by the factor to negative powers. Furthermore, using the 
explicit value for ctq (^2]), the fluctuations become massless at the end of the spectrum, 
indicating that the procedure fails for the energies close to the end of the spectrum. 

In the standard supersymmetry method, it is assumed that uj' ~ A^~^ and only the 
highest order terms in the exponent of the partition function p6| ) are kept. Further, 
all terms proportional to uo in the expressions for the saddle point (^) are neglected. 
As a consequence, the massive modes decouple and are integrated out. The expression 
for the averaged partition function obtained by the above approximations deflnes the 
non-linear zero- dimensional a model. The leading order term of the connected part 
of the two point correlation functions can be calculated from such partition function 



[^1^. If the corrections to that result are needed, or the same correlator have to be 
calculated for cu' ~ 1, other terms have to be considered also, as in ( ^Bj ) and (^). It 
is the interaction between the modes parameterizing the matrix Tq and the massive 
modes, 5P that give rise to the non-leading order terms. 

4 The Correlation Functions 

The correlation functions that have to be calculated are obtained from the second 
derivative of the partition function with respect to the sources, as can be seen from 
(|^ and (|T^). Differentiating and setting the sources to zero, the expression 



whose evaluation is the main aim in this section is obtained: 

I = A{G{E,)G{E,)) = (^(^)) Uo = / VaKuK.^e''^'', (47) 
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where: 



Kpp = —tr^ 



pp 



ip = l,2), A;""' = 



(48) 



B + F = -2Nuj+'tTg To\cxo + N-^/^6P)ToL + 2N^/^uj+'tTg L6P - 

^trg [l - (ao + 200+' Lf] i6P)'. (49) 

The measure is given by (0), B is the part of the exponent that does not contain 
any Grassman variables and F is the Grassman dependent part of the exponent. 
The calculation of the above expression begins with the evaluation of the exponent: 



B 



-2iVcu+'(aoi -ao2)(Ao- Ai) - 

2iVi/2^+' [(Ao - l)(ai - aa) - ^(Ai - - 62)] 



1 r 
2 



1 r 
2 



1 - ((702 - 2u;+'f 



(ai + hi 



(50) 



1 - ((701 + 2c^+ )2J {al + bl) 

F = 2N'/'uj+'{Xo - Ai) [(ai - zhM - Ci){v ' Ci) + («2 - ^b2){p* - C^P - C2)] • 

(51) 

It is understood that Uop take value: aQ^^ depending on which one of the two two-point 
correlators, {G+ {Ei)^^ {E2)) , is calculated. 

Next, the Grassmanian part of the exponent is expanded, it is multiplied by the 
rest of the integrand of (0) and the product expanded in terms of the Grassman 
variables. The integration over the Grassman variables is done using Wegner's the- 
40| . The special case of the theorem needed here is presented in appendix p. 
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The important steps of the integration over the Grassman variables are summarized 
in appendix |C|, leading to the conclusion that among all terms in the expansion of 
the integrand of (0) in terms of the powers of Grassman variables, only two con- 
tribute to the connected correlator, {G{Ei)G{E2)Y: the term with factor ri*r]p*p, 
denoted by: (G'(£'i)G'(£'2))(^=o), and the term with factor r]*r]p* PC1C1C2C2, denoted 
by: {G{E,)G{E,))l^^^y. 

{G{E,)G{E2)Y = (G(Ei)G(E2))^c=o) + {G{E,)G{E2))l^^^y (52) 
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If the appropriate terms of (|47|) are extracted, the integration over the Grassman 
variables can be completed. Using the explicit expression for the measure, (|37[), and 
integrating over 0o and 0i, the following expression for the connected correlator is 
obtained: 



1 roo rl 



2iVw'(Ao-Ai)(cr„i-o-02) 



(53) 



N-''{GiE^)G{E2)) 



c 

(max) 



-1 , , daida2dbidb2liP)Imaxe^ 

4A2(27r)2A V-i V„oo (Ao-Ai)2(ai-26i)2(a2-z62)2 

where /(-P) and i? are given in ( pUf ) and ([50|), respectively, and: 

/maa^ = {ai - ibi){a2 - ib2){Xo - Xi)^ {2N~^ - Au^' {ai + aa - ibi - ^62)^+ 



(54) 



C"oi ~ cro2 



^01-^02 + ^Ar''/'(&l-&2) 



8A^cj+' AoAi 



O"01 



(To2 + iV ^/^(ai-aa) 



fX01-CT02+«iV ^/2(6i-62) 



(55) 



^CJ+ An 



2 

fToi - ao2 + A^"^/^(ai - aa) + 4A^cj+' A, 



O"01 ~ O"02 



The final step is the integration of (^) and (^) over the ordinary variables. The 
crucial problem at this point is the convergence of the integrals. Clearly, from (^OD, 
the integrals over ai, 02, b\ and 62 are Gaussian, therefore convergent. Furthermore, 
since Aq — Ai > over the entire integration range, the integrals in ( |5^ ) and (^^ are 
convergent as long as the real part of — 2A^ct;'((Toi — cro2)(Ao — Ai) is negative. The 
other terms in the exponent ( pOD are small, since they contain the massive fluctuations: 
ai,a2,&i,&2 < A^^^^. Only in that case does the perturbative expansion around the 
saddle point make sense. The convergence properties, as well as the final calculations, 
are given separately for the two regimes in the next two sub-sections. 



4.1 Large u regime 

In the large uj regime (lo' of the order of, but smaller than 1) one looks for the smoothed 
correlators - all terms oscillating with frequencies of the order iV are neglected. (More 
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precisely, the correlator is averaged over the scale few times bigger than A^~^.) The 
idea is that the oscillations are due to the discrete energy spectrum when is finite 
(isolated singularities of the Green's function), which should disappear when N oo 
(the singularities merge to form a branch cut). The integrations over Aq and Ai in each 
of the expressions for the two parts of the connected correlator (|53D and (Q) give four 
terms. The only non-oscillating term comes from the boundary values Aq = 1 = Ai. 
Nevertheless, even the neglected oscillating terms have to be finite in order for the 
procedure to be valid. In order to motivate the forthcoming discussion, one can easily 
see that the connected two point correlator (|52D, remains the same even if some of the 
massive modes, ai, 02, &i or 62, are scaled by an arbitrary complex constant c, provided 
that the real part of the quadratic term in the exponent remains negative. It may 
be surprising, although it is easy to check, that, considering the same equations and 
bearing in mind the relation between fipfip and Ap (|3TD, if both fiofio and fiifii are 
scaled by c, or equivalently Ap 1-^ cAp + 1 — c, the term that survives smoothing, after 
the integrations over Aq and Ai are done, does not explicitly depend on c. Therefore, c 
can be chosen such that Re cc<j+((To^^ is positive infinitesimal so that all terms 

are convergent, while the non- vanishing term is unaffected by the scaling. In effect, 
what has been done is to rotate the contour of the integration of Aq and Ai in the 
complex plane, while keeping the branch-cut structure of the Green's function intact. 
The rotation is justified because the integrands of ( |53D and (|5^ ) are polynomials in 
Aq and Ai, so the only poles are at infinity. Furthermore, from the discussion given 
in section ^ the upper limit of Aq has to be understood as large but finite. It is set 
to infinity only after the limit ^ 00 is taken. 

The above procedure may resemble the one used in the orthogonal polynomials 
method, in which the same connected correlator was calculated at infinity in the 
complex plane and then analytically continued to the real axis, making sure that the 
appropriate branch-cut is obtained [|T3|. Since the correlator in that case is calculated 
at infinity, the oscillatory parts do not exist: the function is smoothed by construction. 
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Similarly, the scaling by an arbitrary constant c given in this article, introduces small 
negative parts in the exponentials of the oscillating terms, setting them to zero. 

The following are the final expressions for the smoothed connected correlator 
{G^ {Ei)G^ {E2)) in the large u) regime, obtained by the straightforward integration 
of (H) and (H). 

N-'{G-'{E,)GHE,)..,,r = ^7,^^ (56) 



4A2 {2Nuj'y 



N-'{G+{E,)G^{E,)) 




(57) 



iV2(a, 



+± 
01 



0", 



02 



) l-( 



01 



+ 2uj 



'\2 



[a, 



02 



\2l2 



5 [1 - (aoi + 2cu')T - 6 [1 - (aoi + 2a;')'] [1 - (cxos - 2a;')'] + [1 - (^02 - 2lu' 



(aoi - ao2)6 [1 - (aoi + 2uj'yf [1 - (^02 - 2uj')^f J 
In order to derive the above relations, the denominator of I{P), (|26D , has been ex- 



panded in terms of the power series. The validity of the expansion is based on the 
fact that only the small quadratic fluctuations are considered, that is, the massive 
modes are smaller than A^^/^ over the entire integration range. 

Using (^ and (j^, the smoothed ensemble-averaged energy density-density cor- 
relation function, in the large N limit and for: Ei — E2 ^ N^, provided the energies 
are not close to the endpoints of the spectrum, is found to be: 

-1 4A2 - E1E2 



C(^l,^2) 



27r2A^2 



{E, - E^fJ^X^ - EfJAX^ - E, 



+ 



40A^ U^jSU -T)-{^f{3SU -T) 
N'^' (^1 - (i?i/2A)yi -\E2/2Xyf 



(58) 



where: 



S = -4(4A'-E,')(4A'-2E2)-4(4A'-£;2)(4A'-2E2)+244(4A'-E2)(4A'-E2') (59) 



T = 2A{2X)\AX- Ef){AX- El)EiE2 
U = AX^~EiE2. 



(60) 
(61) 
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It should be noted that the usual supersymmetry method (zero-dimensional a 
model) does not give a result in this regime. The leading order of the above result is 



identical to the result obtained by Brezin and Zee [0] giving another confirmation of 
the correctness of the above stated convergence prescription. As for the sub-leading 
term, it is questionable if any of the neglected parts, such as higher than quadratic 
fluctuations, may give contribution of the same order. Nevertheless, it is correct 
within the current approximation, namely small quadratic fluctuations around the 
saddle point. 

4.2 Small cj regime 

Let us consider the convergence of {G^{Ei)G^ {E2) first. Again, the convergence is 
achieved if the real part of 2A^ti;'(croi — cro2)(Ao — Ai) is positive. From (^), for small 
positive e and for small but finite u' (of the order A^~^), it follows that: 



(62) 



Therefore the convergence can be achieved by taking the N 00 limit first, so 
that — > while e is still positive, making the above expression positive. Only 
then can e be set to zero giving the required correlator. Using the procedure just 
stated, the terms coming from the upper bound of the integration over Aq in (EBp 



and (|5^) vanish after the limit — 00 is taken. The integration over Ai gives 
the non-exponential terms (from Ai = 1) and the purely oscillatory terms (Ai = 1), 
rendering the integration convergent. Note that in the usual supersymmetry method 



(cr-model) the terms proportional to tu in the expressions for ctq, (|42|) , are neglected 
and consequently the integrals are immediately convergent. The above analysis was 
necessary for the calculations beyond the cr-model, but could also be thought of as a 
clarification of the usual procedure. 

The integrations needed in ( |53D and (Q) are now straightforward, with the fol- 



22 



lowing result: 



4A2 {2NLu'y 



(63) 



iV-2(G+(E0G"(i?2))^™..) 
-1 1 



-167r2 



4A2 (27r)2 1 jv2( 



1 - (aoV + 2cu')2 
-167r2 



02 



- 2cj 



'02 



2cj 



/^2 



(64) 



iTTUJ 



/\2 



0-, 



01 



+ 2uj 



/\2 



(J, 



02 



2cu 



/^2 



,8i7VwVl--B'2 



In order to derive the above relations, the denominator of I{P), (^61) , as well as 
the factor exp{—4:iN^^'^uj' (bi — 62)} appearing after the Ai integration in (|5^ , have 
been expanded in terms of the power series. The validity of the expansions is based 
on the fact that only the small quadratic fluctuations are considered, that is, the 
massive modes are smaller than A^^/2 over the entire integration range. 

Using the expressions for the saddle points (H3) the final expression is obtained: 



iV-2(G+(Ei)G~(E2))^„..) 

-1 



16A2Ar2(l-E'^^ 



E' 



E'- 



e^^^^'v^^^ (65) 



The equations (^) and (^) together constitute the final result for the connected 
two point correlation function {G^ {Ei)G~ {E2)) in the small u regime. While the 
former is a well known result obtained using the standard supersymmetry method, 
among other methods, the latter is the leading order correction to the result. 



5 Conclusion 

One of the methods used in RMT is the supersymmetry method. The usual procedure 
includes some approximations and leads to the zero dimensional cr-model. In this 
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article, the Gaussian unitary ensemble is considered and the supersymmetry method 
is extended to keep some of the previously neglected terms. The partition function 
for the extended method is given in (|i6|). It contains the cx-model terms, as well as 
the lower order terms (quadratic fluctuations around it). 

In the large energy difference regime, the extended supersymmetry method is 
used to calculate the connected energy density- density correlation function (|^. The 
usual degeneracy of the saddle point is hfted fully. Using a symmetry of the model 
(invariance under rescaling of the variables fipfip) the potentially divergent integrals 
can be made finite. As a consequence the smoothed connected energy density- density 
correlation function is calculated, (|58|) . Both the leading order term, which agrees 



with the known result obtained by the method of the orthogonal polynomials, 
and the highest order corrections in the expansion are contained in the result. 
Furthermore, both contributions can be seen to come from terms that are neglected 
in the usual supersymmetry approach. 

As for the small energy difference regime, the degeneracy of the saddle point is only 
slight. The contributions from the (non-degenerate) saddle point, as well as from the 
quadratic fluctuations are taken into account. Convergence can be achieved provided 
that the right limiting procedure is used. As the result, the leading correction to the 



well known result for {G^ {Ei)G~ (£2))'^ is calculated (|65D . 

A natural extension of the work presented here would be to consider other ensem- 
bles (orthogonal, symplectic, . . .). Also, the application of the modified supersymme- 
try method to non-Gaussian ensembles seems worth investigating. 
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A Grassman algebra: definitions and conventions 

Grassman variables, rji and 772 are anti-commuting: 



V1V2 = -V2Vi- (66) 
The complex-conjugate of a Grassman variable is defined as: 

(^^2)* = V*iV*2, V** = -V- (67) 

A general super-vector ip is defined as: 

ip = {zi, r]i, Z2, r]2, ■ ■ •) (68) 

where Zp and rjp are vectors with ordinary and Grassman elements respectively. The 
supervector (p used in this article ((|14D and below) has the form: 

<^ = {<^i, <^l, <^l, <^l), (69) 

where is a A^-dimensional vector consisting of ordinary [a = 0), or Grassman 
(a = 1) elements 'Pp^jj) with yU = (1, . . . , A^). Also, (p is defined as = (p^^ s where 
s = {In, In, -Ijv, Ijv), or explicitly: 

Tp=i^T,^r,-^2\^2)- (70) 

If a and b are ordinary (Bosonic) and a and (3 are Grassman (Fermionic) variables 
or matrices, the super-transpose and Hermitian adjoint of a super-matrix A are: 

a, b, a and (3 are called boson-boson, fermion-fermion, boson-fermion and fermion- 
boson blocks. The super-matrix A was given in what it is called a boson-fermion 
notation. A super-matrix in a general form can be obtained from a super-matrix in the 
boson-fermion block notation by interchanging some of its rows and columns, so that 
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it appears to consist of blocks, each block having the boson-fermion block form. The 
definitions for the super-transpose and Hermitian adjoint are to be applied by putting 
the matrix in boson-fermion notation first and then using the above definitions. 

Elements, or blocks, of both super- vectors and super-matrices used in this article 
have their superscripts denoting grading (0 for ordinary and 1 for Grassman variables), 
while the subscripts denote blocks, each having the boson-fermion block form. 



B Grassman integration and Wegner's theorem 

In this appendix, the integration over the Grassman variables is defined, the Gaussian 
super-integral is given and a special case of Wegner's integral theorem, needed in this 
article, is presented. 

Grassmanian integration is defined as: 

J d7] = 0, J d7]7] = (27r)-^/2 (72) 

It is easy to extend Gaussian integration to include both ordinary and Grassman 
variables: 

/ n idzkdzl / n drjkdvile" (j)^ Acf) = detg'M, (73) 
where = (2;, rf) is a 2A^-dimensional supervector, A is a super-matrix in the boson- 



fermion notation and detg is the graded (super) determinant |^ . The above integral 
is commonly referred to as the Gaussian superintegral. 

A property of the Grassman integration, called Wegner's integral theorem, de- 
serves clarification. Although it is not often explained, it is important for the calcu- 
lations in this article. The motivation for the theorem comes from considering the 
following simple Gaussian superintegral: 

I = j dpdqdr]dr]*e-^^P'+'''^+'^''* = 1 (74) 

On the other hand, performing the following change of variables: 

p = u + CCiu-iv), 7] = Ciu-iv), 
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iq = iv — (*({u — iv), r]* = ({u — iv) (75) 
and including the Jacobian factor: {u — the integral becomes: 

' dudvdCdCiu - ivy^e-^^'^'-^^'l (76) 



Although there is no Grassman dependence in the integrand, the integral is not zero 
because the integration over the ordinary variables is singular. The correct result 
was already given in (|7^. In many cases, when an integral of the form of ( [75| ) is 



encountered, it is neither easy, nor convenient to transform it to the form of ([74|). 
The theorem gives the prescription on how to evaluate such an integral. Instead of 
the full theorem, only a special case of the theorem, relevant for the calculations 
presented in this article, is now presented. A more precise formulation is given in 

Consider an integral of the form: 

/ = Jvix{Q)f{Q), (77) 

where Q is a super-matrix that can be parametrized as a coset space: Q = U^^AU. U 
contains Grassman and A contains ordinary variables, while T>fi{Q) has a singularity 
at the origin of the coset space. / can be expanded as f{Q) = fo + fmax + father, 
where /o and fmax are the terms in the expansion of f{Q) in the powers of Grassman 
variables that contain no Grassman variables and the product of all other terms in 
the expansion. If father goes to zero at the origin, the integral / can be evaluated to 
give: 



1=1 Vn{Q)fmax + /O 



(7J 

A=Ao 



where Aq is the origin of the coset space. 



C Integration Over Grassman Variables 

The integration over Grassman variables in the integral (^^ is outlined in this ap- 
pendix. As a result, the terms in the Grassman expansion of the integrand relevant 
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for the connected correlators, are identified. 

The integration over the Grassman variables of the matrix To {t], t]*, p and p*) 
can be done by expanding the integrand in terms of those variables using the obvious 
notation: 

I = J Po-(/o + fmaxV*'nP*P + father), (79) 

father vauishes at the origin of the Tq coset Aq = 1 = Ai (the 0o,i are not important 
because they decouple). The Wegner's theorem applies, giving: 



Vaf^ax. (80) 

To=l J 

The first term of the above equation can be easily computed using (0), (|^) and 
r 4N'^ 

/ VPI{P)fo = -^(^01 + 2u;+ ){ao2 - 2uj+ ). (81) 

J To=i y 

On the other hand, {G{Ep)) can be calculated starting from (|12D and using (H), (^8]) 
and (^9]). Wegner's theorem is used again to show that the only contribution comes 
from the term that has no Grassman variables and gives: 

{G{E,)) = -j{a,,-{-y2u;+'). (82) 



Therefore, the first term in (^0]) (divided by 4) is exactly the disconnected part of 
the correlator {G{Ei)){G{E2)). Therefore the connected part (^), is simply: 

{G{E,)G{E2)r = \j VaUa. ^ (83) 

Next, P is integrated over the Grassman variables contained in 5P: Ci, C* > C25 and 
Q- The procedure is the same as above, only this time it has to be done separately for 
5Pi and 5P2 because they are independent and both of them have a coset structure. 
The origins of the respective coset spaces are given by: ai = = 61 and 02 = = 62, 
so that P can be expanded as: 

= jVa [g,^ + g,,,^^C;Cp + Qvother] ' (84) 
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in obvious notation, for both p = 1 and p = 2. Both gp^^f^^^ vanish at the respective 
origins, so Wegner's theorem apphes. Consequently, the connected correlation func- 
tion receives only two contributions: one from the term with factor 7]*rip* p evaluated 
at the origin of 5P coset (denoted by subscript (C = 0) ) and the other from the 
term with all Grassman variables (denoted by subscript (max)). Therefore the above 
statement can be written as: 

{G{E,)G{E,)Y = (G(Ei)G(E2))^c=o) + {G{E,)G{E,))l^^^y (85) 

The above equation specifies how to perform the integration over the Grassman 
variables needed in section |[ 
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